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O ■ ABSTRACT 



Renormalization group, and in particular its Quantum Field Theory imple- 
mentation has provided us with essential tools for the description of the phase 
transitions and critical phenomena beyond mean field theory. We therefore review 
the methods, based on renormalized quantum field theory and renormaliza- 
tion group, which have led to a precise determination of critical exponents of 
the N- vector model and of the equation of state of the 3D Ising model ||. 
These results are among the most precise available probing field theory in a non- 
perturbative regime. 

^ ■ Precise calculations first require enough terms of the perturbative expansion. 

However perturbation series are known to be divergent. The divergence has been 
characterized by relating it to instanton contributions. The information about 
large order behaviour of perturbation series has then allowed to develop efficient 
"summation" techniques, based on Borel transformation and conformal mapping 

We first discuss exponents and describe our recent results []2[]. Compared to 
exponents, the determination of the scaling equation of state of the 3D Ising 
model involves a few additional (non-trivial) technical steps, like the use of the 
parametric representation, and the order dependent mapping method. From the 
knowledge of the equation of state a number of ratio of critical amplitudes can 
also be derived. 

Finally we emphasize that few physical quantities which are predicted by renor- 
malization group to be universal have been determined precisely, and much work 
remains to be done. Considering the steady increase in the available computer 
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resources, many new calculations will become feasible. In addition to the infi- 
nite volume quantities, finite size universal quantities would also be of interest, 
to provide a more direct contact with numerical simulations. Let us also mention 
dynamical observables, a largely unexplored territory. 



1 Introduction: the effective 4 field theory 

Second order phase transitions are continuous phase transitions where the cor- 
relation length diverges. Renormalization group (RG) arguments ||, as well as 
an analysis, near dimension four, of the most IR divergent terms appearing in 
the expansion around mean field theory || , indicate that such transitions possess 
universal properties, i.e. properties independent to a large extent from the details 
of the microscopic dynamics. Moreover all universal quantities can be calculated 
from renormalizable quantum field theories which are local when the interactions 
are short range. For an important class of physical systems and models one is 
led to a (/> 4 -like euclidean field theory with O(N) symmetry. Among those let 
us mention statistical properties of polymers, liquid-vapour, binary mixtures, 
superfluid Helium, ferromagnetic transitions... We explain here how critical ex- 
ponents and other universal quantities have been calculated with the help of 
renormalization group ideas and their implementation within the quantum field 
theory context 0. 

The effective quantum field theory. Wilson's renormalization group shows that, 
at least in the neighbourhood of dimension four, the universal properties of statis- 
tical models which short range interactions and O(N) symmetry can generically 
be described by an effective euclidean field theory, with an action Ti (</>) 

H {(!>)= J d d x{i[V^(x)] 2 + ^ 2 (x) + ^ A 4 - d (^ 2 (x))) 2 }, (1.1) 

where r is a regular function of the temperature T. The parameter A has the 
dimension of a mass and corresponds to the inverse microscopic scale. It also 
appears as a cut-off in the Feynman diagrams of the perturbative expansion. 

For some value r c = r(T c ) the correlation length £ diverges (the physical mass 
m = l/£ vanishes). Near the critical temperature T c , in the critical domain 
\T — T c \ oc \r — r c \ <C A 2 , | (0) | <C A^ -2 )/ 2 , the mass remains small, ro<C A. In 
this limit a universal behaviour of thermodynamic quantities is expected. The 
study of the critical domain thus reduces to the study of the large cut-off be- 
haviour, i.e. to renormalization theory and the corresponding renormalization 
group. However, in the traditional presentation of quantum field theory in the 
context of particle physics, the dependence of the parameters of the action as a 
function of the cut-off A is determined by the condition that renormalized cor- 
relation functions should have a finite cut-off limit. Here instead the coefficient 
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goA 4-d /4! of the (/> 4 interaction has a dependence on A given a priori. In partic- 
ular in the dimensions of interest, d < 4, the "bare" coupling constant diverges 
for A — > oo, a reflection of the IR instability of the gaussian fixed point, though 
the field theory, being super-renormalizable, requires only a mass renormaliza- 
tion. 

The divergence of the bare coupling constant. One solution to the problem of 
the large coupling constant is provided by Wilson-Fisher's famous e-expansion 
||. One defines, at least in perturbation theory, the field theory in arbitrary 
complex space dimension d. Setting d = 4 — s one then expands both in go 
and e. The large cut-off divergences become logarithmic as in four dimensions, 
and can be removed by standard renormalizations. Renormalization group (RG) 
equations for correlation functions follow, from which scaling laws can be derived, 
and which lead to calculations of physical quantities as power series in e. 

The most convenient RG equations are homogeneous differential equations, 
first derived for the critical (massless) theory ||, and then generalized to the 
whole critical domain r ^ r c ||10|1 . They correspond to a scheme in which the 
massless (or critical) theory has a perturbative expansion in go and therefore rely 
on the double series expansion in go and e. An alternative formalism is based 
on the Callan-Symanzik (CS) (inhomogeneous) equations for the massive theory 
IT]. One introduces correlation functions for a renormalized field 0r = Z -1 / 2 0, 



expressed in terms of the physical mass m or inverse correlation length £ _1 = m 
and a renormalized coupling constant g. They are implicitly defined by the 
conditions (after factorization of the obvious group indices): 



T£>(p;m,g) = m 2 + p 2 + O (p 4 ) , (1.2a) 

R 



f ( £ ) (p l = 0;m 1 g)=m A - d g 1 (1.2b) 



where the functions are the Fourier transforms of so-called (1PI) correlation 

functions , the coefficients of the expansion of the thermodynamic potential 
in terms of the local field expectation value (4>n{x)). 

To prove scaling laws from the CS equations it remains necessary to expand 
in g and e. The form of the CS /3-function in dimension d is 

(3(g) = -{A-d)g + a 2 (d)g 2 + ■ ■ ■ , a 2 (d) > . 

For e = 4 — d small it has a zero g* of order e and one recovers the principle of 
the £-expansion. 

However in this framework, because the theory now is massive, the perturba- 
tive expansion is IR finite in any dimension and the CS equations are always 
satisfied. Moreover one generally assumes that the results established for the 
e-expansion remain valid at finite e. Therefore Parisi |T3] has suggested work- 
ing at fixed dimension d < 4, in the massive theory (the massless theory is IR 
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divergent). In contrast with the e-expansion however, at fixed dimensions three 
or two no small parameter is available. Therefore a precise determination of 
g* and all other physical quantities depends on the analytic properties of the 
series, in addition to the number of terms available. A semi-classical analy- 
sis, based on instanton calculus, unfortunately indicates that perturbation in 
quantum field theory is always divergent. Therefore to extract any information 
from perturbation theory a summation method is required. Note finally that at 
fixed dimension, at any finite order, universal quantities are given by expressions 
which depend on the renormalization scheme, in contrast with the results of the 
£-expansion. 

This approach, which is more cumbersome, has eventually been tried for prac- 
tical reasons: it is easier to calculate Feynman diagrams in dimension three than 
in generic dimensions, and thus more perturbative orders could be obtained. 

While this approach is a natural extension of the ^-expansion, its interpretation 
directly at fixed dimension is worth discussing. At fixed dimension d < 4, in the 
critical domain (condition which implies shifting r by r c and thus performing a 
mass renormalization) , all terms in the perturbative expansion have a finite large 
cut-off limit at uq = goA 4 ~ d fixed because the theory is super-renormalizable. 
This means physically that the initial parameters of the field theory are first 
tuned to remain artificially close to the unstable uq = gaussian fixed point. 
Indeed the true expansion parameter is a dimensionless quantity and therefore 
is proportional to goA 4 ~ d /m 4: ~ d . To keep this parameter constant when m the 
physical mass goes to zero for T — > T c (r — ► r c ) in the high temperature phase, 
one must vary the dimensionless parameter go as 

gl /{A ~ d) oc m/A oc 1/fA. 

The relevant theory, however, corresponds to the infinite uq limit. One is then 
confronted with a serious technical problem: perturbation theory is finite in the 
critical domain but one is interested in the infinite coupling limit, where obviously 
the perturbative expansion is no longer useful. 

One thus introduces the renormalized coupling constant g and the field renor- 
malization defined by the conditions (|1.2j) . Then 

u = g A^ = m^G(g), (3(g) = (d- A)G{g)/G'{g). (1.3) 

When the initial coupling constant uq becomes large the new dimensionless cou- 
pling g has a finite limit provided the /3-function has an IR stable zero g*\ 

Uo _, oo => /3(g*) =0 andwE p'(g*) > . 

In such a situation the renormalized coupling g is a more suitable expansion 
parameter then uq. 
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The relations (\1.3j ) then imply that at go fixed, A/ra — >• oo, 

<7-<7*cx(m/A) w . (1.4) 

To the field renormalization Z(g) is associated the RG function r](g) 

rj(g) = /3(g)Z'(g)/Z(g). (1.5) 

Integration of the equation ( |1.5|) yields the behaviour of Z(g) near g* 

Z{g)<x(g*-gY'" a (to/A)", (1.6) 

where the exponent = ?7(<7*) characterizes the field anomalous dimension d^: 
2d(f) = d — 2 + r\. This singular behaviour, consistent with the scaling proper- 
ties derived from the ^-expansion, explains why even in dimension d < 4 the 
introduction of a field renormalization is necessary. 

The approach has sometimes be questioned, because it involves this double 
limit, but the final results and their comparison with other data have shown the 
consistency of the method. 

A last remark: in this framework the mass parameter m has still to be related 
to the temperature. RG equations show that it is singular at T c and behaves for 
t = r — r c — > 0_|_ as 

m oc t v oc (T - T c f , (1.7) 
where v is the correlation length exponent. 

2 Critical exponents and series summation 

Critical exponents are the most studied quantities in the theory of phase tran- 
sitions, because they characterize universality classes and are easier to calcu- 
late. They have been extensively used to test RG predictions by comparing 
them with other results (experiments, high or low temperature series expansion, 
Monte-Carlo simulations) ]13[] . The first precise determination of the exponents 
of the O(N) symmetric A^-vector model has been reported by Le Guillou-Zinn- 
Justin UJ, using the six-loop series for RG functions calculated by Nickel et al 



14| . The summation methods used in []T|] were based on a Borel transformation 
and conformal map. The same ideas have later been applied to the e-expansion 
15| when five loop series became available |l(|, and more recently to the equa- 



tion of state by Guida-Zinn- Justin [Q. With time the method has been refined 
and the efficiency improved by various tricks but the basic principles have not 
changed. 
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2.1 Borel transformation and conformal mapping 
Let F(g) be any quantity given by a perturbative series 



F{g) = Y J F k9 k . (2.1) 

k=0 

A semi-classical analysis of barrier penetration effects for negative coupling (in- 
stantons) teaches us that in the 0§ field theory for large order k the coefficients 
Fk behave like 

F k oc k s {-afk\. 

k — >00 

The value of a > has been determined numerically (while s is known analyti- 
cally). One then introduces B(b, g), the Borel (rather Borel-Leroy) transform of 
F(g), which is defined by 

^») = E^i/. < 2 - 2 > 

where b is a free real parameter [b > — 1). Formally, i.e. in the sense of series 
expansion, F(g) can be recovered from 

POO 

F(g)= / t 6 e-*fl(&,c#)d*. (2.3) 
Jo 

Using the large order estimate in (|2.2| ) one verifies that B(b, z) is analytic at 
least in a circle of radius 1/a and its singularity closest to the origin is located 
at z = —1/a. Therefore B(b,z), in contrast with F(z), is determined by its 
series expansion. However, the relation (|2.3|) can be extended from a relation 
between formal series to a relation between functions only if B(b, z) is analytic in 
a neighbourhood of the real positive axis. In the case of the 0| field theory such 
a property has been proven in constructive field theory (as well as the property 
that the function F(g) is indeed by given by ( |2.3|) ) |T7| . Moreover it is necessary 
to construct an analytic continuation of the function B(b, z) from the circle to 
the real positive axis. Considerations about more general instanton contributions 
strongly suggest that the Borel transform actually is analytic in a cut-plane, the 
cut being on the real negative axis, at the left of —1/a. Therefore an analytic 
continuation is provided by a conformal map of the cut-plane onto a circle: 

. , y/1 + az - 1 
z h-> u(z) = —=^=-— . (2.4) 
Vl + za + 1 

The function B[b, z(u)\ is then given by a series in powers of u convergent in 
the cut-plane. The corresponding (hopefully convergent) series expansion for 
function F{g) takes the form 

oo /»oo 

F(jg) = TTu k {b) / t b e"' [u(gt)] k dt . (2.5) 
k=o Jo 
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The parameter 6, as well as a few other parameters introduced in variants, are 
used to improve the apparent convergence and test the sensitivity of results 
to their variations. Moreover the value of b has to stay within a reasonable 
range around the value s predicted by the large order behaviour. Finally the 
summation method is expected to be efficient mainly when the available co- 
efficients Ffc behave already as predicted by the asymptotic large order esti- 
mate. 

2.2 Exponents 

The values of critical exponents obtained from field theory have remained after 
about twenty years among the most precise determinations. Only recently have 
consistent, but significantly more precise, experimental results been reported in 
low gravity superfluid experiments [[IT] . Also the precision of results coming from 
high temperature expansions [|18| -p9f and various numerical simulations |30|-f40 



on the lattice has kept steadily improving. 

Recently seven-loop terms have been obtained for < N < 3 for two of the 
three RG functions, related to the <p and (p 2 dimensions, by Murray and Nickel 



42fl . These terms, together with some improvement in the summation methods, 
have led to the new slightly more precise values of g* and critical exponents 
displayed in table | (gfc = (N + 8)g* /48tt) (Guida-Zinn- Justin @). Among the 
exponents 7, v, rj, (3, a, only two are independent, for example 

7 = 1/(2-7/), /?= 1^(1 + 7/), a = 2-3u, 

but they are calculated independently to test the precision of the summation 
method. 

The main improvements concern the exponent rj which was poorly determined, 
and the lower value of 7 for A = (polymers) . 

In the framework of the ^-expansion, since the series used earlier were affected 
by a small error at order e 5 , the values have also been updated (table |2|). Two 
kinds of results are reported, free means simple summation as above, be means 
that the known values in two dimensions have been incorporated in the summa- 
tion procedure. It is gratifying that the overall consistency between the two set 
of values coming from 3D series and ^-expansion has improved. 

Finally 0(4) results of interest for simulations of the Higgs phase transition at 
finite temperature, obtained from six loop series, have been added (table |3|) . 

3 The scaling equation of state 

Let us first recall a few properties of the equation of the state in the critical 
domain, in the specific case A = 1 (Ising-like systems), at d < 4. 

The equation of state is the relation between magnetic field H, magnetization 
M = (</>) (the "bare" field expectation value) and the temperature which is 
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Table @ 

Critical exponents of the 0(N) models from d = 3 expansion J^/. 



N 





1 


2 


3 


~ * 


1.413 ±0.006 


1.411 ±0.004 


1.403 ±0.003 


1.390 ±0.004 


<?* 


26.63 ±0.11 


23.64 ±0.07 


21.16 ± 0.05 


19.06 ±0.05 


7 


1.1596 ±0.0020 


1.2396 ±0.0013 


1.3169 ±0.0020 


1.3895 ±0.0050 


7 / 

Is 


0.5882 ±0.0011 


0.6304 ±0.0013 


0.6703 ±0.0015 


0.7073 ± 0.0035 


V 


0.0284 ±0.0025 


0.0335 ± 0.0025 


0.0354 ±0.0025 


0.0355 ±0.0025 




0.3024 ± 0.0008 


0.3258 ±0.0014 


0.3470 ±0.0016 


0.3662 ±0.0025 


a 


0.235 ±0.003 


0.109 ±0.004 


-0.011 ±0.004 


-0.122 ±0.010 




0.812 ±0.016 


0.799 ±0.011 


0.789 ±0.011 


0.782 ±0.0013 


6 = ujv 


0.478 ±0.010 


0.504 ± 0.008 


0.529 ±0.009 


0.553 ±0.012 






Table g 








Critical exponents of the O(N) models from e-expansion ||/. 


N 





1 


2 


3 


7 (free) 
7 (be) 


1.1575 ±0.0060 
1.1571 ±0.0030 


1.2355 ±0.0050 
1.2380 ±0.0050 


1.3110 ±0.0070 
1.317 


1.3820 ±0.0090 
1.392 


v (free) 
(be) 


0.5875 ±0.0025 
0.5878 ±0.0011 


0.6290 ±0.0025 
0.6305 ± 0.0025 


0.6680 ± 0.0035 
0.671 


0.7045 ± 0.0055 
0.708 


r\ (free) 
?7 (be) 


0.0300 ±0.0050 
0.0315 ±0.0035 


0.0360 ± 0.0050 
0.0365 ± 0.0050 


0.0380 ± 0.0050 
0.0370 


0.0375 ±0.0045 
0.0355 


/3 (free) 
P (be) 


0.3025 ±0.0025 
0.3032 ±0.0014 


0.3257 ±0.0025 
0.3265 ±0.0015 


0.3465 ± 0.0035 


0.3655 ±0.0035 




0.828 ±0.023 


0.814 ±0.018 


0.802 ±0.018 


0.794 ±0.018 





0.486 ±0.016 


0.512 ±0.013 


0.536 ±0.015 


0.559 ±0.017 



represented by the parameter t = r — r c oc T — T c . It is related to the free energy 
per unit volume, in field theory language the generating functional T(<j)) of 1PI 
correlation functions restricted to constant fields, i.e the effective potential V, 
V(M) = r(M)/vol., by H = dV/dM. In the critical domain the equation of 
state has Widom's scaling form (5 = (d + 2 — rj)/(d — 2 ± rj)) 

H(M, t) = M J /(t/M^), (3.1) 

a form initially conjectured and which renormalization group has justified. 
One property of the function H (M, t) which plays an essential role in the 
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Table | 

Critical exponents in the 0(4) models from d = 3 and e-expansion j^j. 





d = 3 


e : free, be 


~ * 
9Ni 


1.377 ±0.005 




g* 


17.30 ±0.06 




7 


1.456 ±0.010 


1.448 ±0.015 , 1.460 




0.741 ±0.006 


0.737 ±0.008 , 0.742 


V 


0.0350 ± 0.0045 


0.036 ±0.004, 0.033 





0.3830 ± 0.0045 


0.3820 ± 0.0025 


a 


-0.223 ±0.018 


-0.211 ±0.024 




0.774 ±0.020 


0.795 ± 0.030 


e 


0.574 ±0.020 


0.586 ±0.028 



analysis is Griffith's analyticity: it is regular at t = for M > fixed, and 
simultaneously it is regular at M = for t > fixed. 

Amplitude ratios. Universal amplitude ratios are numbers characterizing the 
behaviour of thermo dynamical quantities near T c . Several of them commonly 
considered in the literature can be directly derived from the scaling equation of 
state. Let us just give two examples. 

The singular part of the specific heat, i.e. the 4> 2 two-point correlation function 
at zero momentum, behaves like 

C H =A ± \t\~ a , tocT-T c ^±0. (3.2) 

The ratio A + /A~ then is universal. 

The magnetic susceptibility x i n zero held, i.e. the <p two-point function at 
zero momentum, diverges like 

x = C ± \t\~ 1 1 t^±0. (3.3) 

The ratio C + /C~ is also universal. 

The e-expansion. The first results concerning the scaling equation of state 
have been obtained within the framework of the e = 4 — d expansion. 

The e-expansion of the scaling equation of state has been determined up to 



order e for the general O(N) model ||43|| , and order e for N = 1 |j44| . We give 



here the function f(x = t/M 1 ^) of eq. (|3.1| ) for iV = 1 up to order e 2 to display 
its structure, 

f(x) = 1 + x + \e{x + 3)L + e 2 [±[x ± 9)L 2 + §^(x + 3)L] ± O (e 3 ) , (3.4) 
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with L = log(x + 3). 

The expression ( |3.4| ) is not valid for x large, i.e. for small magnetization M. 
In this regime the magnetic field H has a regular expansion in odd powers of M. 
It is thus convenient to express the equation of state in terms of another scaling 
variable z oc x"' 3 because 

z oc Mt~P. (3.5) 
The equation of state then takes the form 

Hoct 05 F(z), (3.6) 

where the relation between exponents 7 = (3(6 — 1) has been used. Substituting 
into eq. ( |3.4j) x = xqz~ 1 ^ (the constant xq takes care of the normalization of z) 
and expanding in e one finds at order e 2 for the function (|3.6| ) 

F(z) = F (z) + eF^z) + e 2 F 2 (z) + 0(e 3 ), 

with 

F = z+±z 3 (3.7a) 
Fx = ±(-z 3 + L(2z + z 3 )) 

F 2 = (~50z 3 + L(100z - Az 3 ) + L 2 (18z + 27z 3 )) (3.76) 



and L = log(l + z 2 /2). 

Within the framework of the formal ^-expansion one can easily pass from one 
expansion to the other. Still a matching problem arises if one wants to use 
the e-expansion to determine the equation of state for d = 3, i.e. e = 1. One 
is thus naturally led to look for a uniform representation valid in both limits. 



Josephson-Schofield parametric representation |45[] has this property. 



4 Parametric representation of the equation of state 

In order to implement both Griffith's analyticity and the scaling relation, one 
parametrizes the equation of state in terms of two new variables R and 6*, setting: 

M = moR^e , 

t = R(l-6 2 ), (4.1) 
H = h R p5 h(9) , 



where ho, mo are two normalization constants. We choose ho such that 



h(9) = e + o(e 3 ) 
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In terms of the scaling variables x of eq. (|3.4|) or z from eq. ( |3.5|) this parametriza- 
tion corresponds to set 

z = p9/(l - 9 2 f, 9>0, (4.2) 

x = x oP - 1 ^(i-e 2 )e- 1 ^, (4.3) 

where p is some other positive constant. 

Then the function h(9) is an odd function of 9 which from Griffith's analyticity 
is regular near 9 = 1, which is x small, and near 9 = which is x large. It vanishes 
for 9 = 9q which corresponds to the coexistence curve H = 0, T < T c . 

Note that the mapping (|4.2| ) is not invertible for values of 9 such that z'(9) = 0. 
The derivative vanishes for 9 = 1/ yjl — 2/3) ~ 1.69. One has to verify that this 
value is reasonably larger than 9q, the largest possible value of 9. 

Finally it is useful for later purpose to write more explicitly the relation be- 
tween the function F(z) of eq. ( |3.6|) and the function h(9): 

h(9) =p~ x (l - 9 2 f 5 F (z{9)) . (4.4) 

Expanding both functions 

F(z)=z + ±z 3 + Y,F2i+iZ 2l+1 , (4.5) 

1=2 

h(9)/9 = l + J2h2i+i9 21 , (4.6) 
i=i 

one finds the relations 

^3 = \ P 2 ~ 7 (4.7a) 

h 5 = | 7 (7 - 1) + |(2/3 - 7 )P 2 + F 5 p 4 (4.76) 
kr = ±7(7 - 1)(7 - 2) + A(2/3 - 7 )(2/3 - 7 + l)p 2 

+ (4/3 - 7 )iV + F 7 p 6 (4.7c) 



From the parametric representation of the equation of state it is then possible to 
derive a representation for the singular part of the free energy per unit volume 
as well as various universal ratios of amplitudes. 

Parametric representation and e-expansion. Up to order e 2 the constant mo 
(or p) can be chosen in such a way that the function h(9) reduces to: 

h(9) = 9 (1 - f 9 2 ) + O (e 2 ) . (4.8) 

The minimal model in which h(9) is approximated by a cubic odd function of 9 
is called the linear parametric model. At order e 2 the linear parametric model 
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is exact, but at order e 3 the introduction of a term proportional to 9 5 becomes 
necessary [4"4",4"6|. One finds: 

h{6) = 6(1 + h 3 6 2 + M 4 ) + O (e A ) , (4.9) 

with 



(4.10) 



where A is the constant 

A = ^'(1/3) - |tt 2 = 1.17195361934 .... (4.11) 
The function h(9) vanishes on the coexistence curve for 6 = 9q: 

Note that h% and thus 9q are determined only up to order e 2 . It follows 
p 2 = 6( 7 + h 3 ) = 2 (1 + \e + j^e 2 ) = 3.13 ± 0.13, 

because h% is determined only up to order e 2 , 

Remark. In the more general 0(N) case, the parametric representation also 
automatically generates equations of state with satisfy the required regularity 
properties, and thus leads to uniform approximations. However for N > 1 the 
function h{6) still has a singularity on the coexistence curve, due to the presence 
of Goldstone modes in the ordered phase and has therefore a more complicated 
form. The nature of this singularity can be obtained from the study of the 
non-linear cr-model. It is not clear whether a simple polynomial approximation 
would be useful. For N = 1 instead, one expects at most an essential singularity 
on the coexistence curve, due to barrier penetration, which is much weaker and 
non-perturbative in the small e- or small (/-expansion. 

5 Perturbative expansion at fixed dimension three 

We now discuss the calculations based on a perturbative expansion at fixed 
dimension d = 3. Five loop series for the renormalized effective potential of the 
03 theory have been first reported in Bagnuls et al. p7j , but the printed tables 



contain some serious misprints. These have been noticed by Halfkann and Dohm 
who have published corrected values [§8j . These five-loop calculations have only 



been performed for N = 1, because they are much more difficult for N ^ 1 due 
to the presence of two lengths, the correlation lengths along the applied field and 
transverse to it. 
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5.1 General remarks 

The general framework again is the massive theory renormalized at zero mo- 
mentum. The correlation functions of the renormalized field 0r = (j)/\fZ 
are fixed by the normalization conditions (|1.2| ). Eventually the renormalized 
coupling constant g has to be set to its IR fixed point value g*. 

The conditions ( |1.2| ) imply that the free energy T expressed in terms of the 
"renormalized" magnetization i.e. the expectation value of the renormalized 
field <p — (0r), has a small (p expansion of the form (in d dimensions) 

F(<p) = ^(0) + imV 2 + i m 4 ~ d gip 4 + O (tp 6 ) . (5.1) 

It is important to remember that the finite field renormalization Z(g) is singular 
at g* (equation (|L6|)). It follows that <p oc 

It is convenient to introduce the rescaled variable z 

z = ipm^-^' 2 ^ oc M/m {d - 2+ri)/2 oc M/t p , (5.2) 

(equation ( p..7|) ) and set 

m 

F( i p)-f(0) = —V(z,g). (5.3) 
9 

Taking into account the definitions of section ^ we conclude that the equation 
of state is related to the derivative F of the reduced effective potential V with 
respect to z 

F{z i9 )=™£&, (5.4) 

by 

H oc t^ s F(z). 

Ising symmetry implies that F is an odd function of z 

F(z,g) =z + l-z 3 + J2F2i+i(g)z 2l+1 . (5.5) 

1=2 



5.2 The problem of the low temperature phase 

To determine the equation of state in the whole physical range, or universal 
ratios of amplitudes, a new problem arises. In this framework it is more diffi- 
cult to calculate physical quantities in the ordered phase because the theory is 
parametrized in terms of the disordered phase correlation length £ = m~ l oc 
(T — T c )~ u which is singular at T c (as well as all correlation functions normalized 
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as in ( |1 . 2p ) . For example at one-loop order for d = 3 the scaling function F(z, g) 
(eq. (|5.4|) ) is given by: 

g) = z + - ^ [(1 + z 2 /2) 1/2 - 1 - z 2 /A\ (5.6) 

= * + J* 3 + 2^9^ - + 0(A 

where the subtractions, due to the mass and coupling normalizations, are deter- 
mined by the conditions (|1.2| ). This expression is adequate for the description 
of the disordered phase. However when t goes to zero for fixed magnetization, 
i.e. m — > at <p fixed, then z —>■ oo as seen in eq. ( |3.5| ). Thus all terms in the 
loopwise expansion become singular. In this limit one knows from eq. ( [3.1| ) that 
the equation of state behaves like 

H(M, t = 0) oc M s F(z) oc z s . (5.7) 

In the framework of the e-expansion the scaling relations (and thus the limiting 
behaviour (|5.7|) ) are exactly satisfied order by order. Moreover the change to the 
variable x oc z~ 1 ^ (more appropriate for the regime t — > 0) gives an expression 
for f(x) oc F(x~^)x^ s that is explicitly regular in x = (Griffith's analyticity) : 
the singular powers of log x induced by the change of variables cancel non trivially 
at each order, leaving only regular corrections. 

The situation changes when one deals with the perturbation theory in d = 
3 dimensions: scaling is not satisfied for generic values of g but only at g*. 
Consequently scaling properties are not satisfied order by order in g. In particular 
the change to the Widom function f(x) will introduce singular powers of x that 
violate Griffith's analyticity. An analogous problem arises if one first sums the 
series at g = g* before changing to the variable x. In this case the singular 
contributions (in the form of powers of x) do not cancel, as a result of unavoidable 
numerical summation errors. 

Several approaches have been proposed to solve the problem of continuation to 
the ordered phase. A rather powerful method, motivated by the results obtained 
within the ^-expansion scheme, is based on the parametric representation . 

Parametric representation and order dependent mapping ( ODM). The problem 
that one faces is the following: to reach the ordered region t < one must cross 
the point z = oo. However we know from Griffith's analyticity that F(z)z~ s 
is regular in the variable z~ x l^ . This property is naturally satisfied within the 
parametric representation. One thus introduces an new variable 9 and an auxil- 
iary function h(9) defined as in ( [4.2| , f4.4| ): the exact function h(9) is regular near 



9 = 1 (i.e. z = oo) and up to the coexistence curve. 

The approximate h{9) that one obtains by summing perturbation theory at 
fixed dimension, instead is still not regular, because the singular terms generated 



by the change of parametrization (|4.2[) at 9 = 1 do not cancel exactly due to 
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summation errors. The last step involves a Taylor expansion of the approximate 
expression of h(9) around 9 = and a truncation of the expansion, to enforce 
regularity. A question then arises, to which order in 9 should one expand? Since 
the coefficients of the 9 expansion are in one to one correspondence with the 
coefficients of the small z expansion of the function F(z, g*), the maximal power 
of 9 in h(9), should be equal to the maximal power of z whose coefficient can 
be determined with reasonable accuracy. Indeed although the small z expansion 
of F(z) at each finite loop order in g contains an infinite number of terms, the 
evaluation of the coefficients of the higher powers of z is increasingly difficult. 
The reasons are twofold: 

(i) The number of terms of the series in g required to get a precise estimate of 
F[ increases with / because the large order behaviour sets in later. 

(ii) At any finite order in g the function F(z) has spurious singularities in the 
complex z plane (see e.g. eq. ( |5.6| ), z 2 = —2) that dominate the behaviour of the 



coefficients F[ for I large. 

In view of these difficulties one has to ensure the fastest possible convergence 
of the small 9 expansion. For this purpose one uses the arbitrary parameter 
p in eq. ( |4.2|) : one determines p by minimizing the last term in the truncated 
small 9 expansion, thus increasing the importance of small powers of 9 which are 
more precisely calculated. This is nothing but the application to this particular 
example of the series summation method based on ODM ||49|| . 

6 Numerical results 

The calculation proceeds in two steps; first one determines the first coefficients 
of the small field expansion, then one introduces the parametric representation. 

6. 1 The small field expansion 

The determination of the coefficients i^z+i of the small z (small field) expan- 
sion of the function F(z) relies on exactly the same method as for exponents, 
i.e. Borel-Leroy transformation and conformal map. In figure [T] the behaviour 
of F5 in terms of the Borel-Leroy parameter b is displayed. Increasing flatness 
of the curves when k increases, i.e. increasing insensitivity to the parameter b, 
supports the hypothesis that the method indeed converges. 

Because the asymptotic regime sets in later when I increases, the efficiency of 
the summation correspondingly decreases, as expected. Table § contains the re- 
sults of || together with other published estimates of the coefficients of the small 
z expansion of F(z) coming from high temperature series [5C ,[52,53], simulations 



54 1 and derivative expansion of the Exact RG |55|,|56 . 
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Fig. 1 The summed coefficient F$ as a function of the Borel-Leroy parameter b for 
successive orders k. 



Table | 
Equation of state. 





2* 




F 7 x 10 4 


F 9 x 10 5 


e-exp., H 
e— exp., 1 
d = 3,g| 
HT [50 
HT [52 
HT [53 
MC [54 
ERG [5t 


§ 
1 
1 


23.3 
23.4 ±0.1 
23.64 ±0.07 
24.45 ±0.15 
23.72 ± 1.49 
23.69 ± .10 
23.3 ±0.5 
20.72 ±0.01 


0.0177 ±0.0010 
0.01715 ±0.00009 
0.01711 ±0.00007 
.017974 ± .00015 
0.0205 ± 0.0052 
.0168 ±0.0012 
0.0227 ±0.0026 
0.01719 ±0.00004 


4.8 ±0.6 

4.9 ±0.6 
4.9 ±0.5 

5.4 ±0.7 
4.9 ±0.1 


-3.3 ±0.3 
-5.5 ±4 
-7 ±5 

-2.3 ± 1.1 
-5.2 ±0.3 



6.2 Parametric representation 
One then determines by the ODM method the parameter p and the function 



h(9), as explained in section |5T2| . One obtains successive approximations in the 
form of polynomials of increasing degree for h(6). At leading order h(9) is a 
polynomial of degree 5, whose coefficients are given by the relations ( fl.7| ): 

h(9) = 9[l + h 3 (p)9 2 + h 5 (p)9 4 }. (6.1) 



For the range of admissible values for F5 the coefficient of 9 5 given by eq. (|4.7fc|) 
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has no real zero in p. It has a minimum instead 

1 



P 2 = Pl = j^r(l-W)- (6.2) 



Substituting this value of p into expression ( pTT|) one obtains the first approxi- 
mation for h(9). At next order one looks for a minimum pi of \h-j(p)\. One finds 
a polynomial either of degree 5 in 9, when hi has a real zero, or of degree 7 when 
it has only a minimum. It is not possible to go beyond hg(p) because already Fg 
is too poorly determined. 

Note that one here has a simple test of the relevance of the ODM method. 
Indeed, once h(9) is determined, assuming the values of the critical exponents 7 
and /?, one can derive the corresponding function F{z). It has an expansion to all 
orders in z. As a result one obtains a prediction for the coefficients -F2Z+1 which 
have not yet been taken into account to determine h{9). The relative difference 
between the predicted values and the ones directly calculated gives an idea about 
the accuracy of the ODM method. Indeed from the values F5 = 0.01711,7 = 
1.2398,/? = 0.3256, one obtains 

F 7 = 4.83 x 10 -4 , F 9 = -3.2 x 10 -5 , F n = 1.4 x 10 -7 

We see that the value for F 7 is quite close to the central value one finds by direct 
series summation, while the value for Fg is within errors. This result gives us 
confidence in the method. It also shows that the value of Fg obtained by direct 
summation contains little new information, it provides only a consistency check. 
Therefore the simplest representation of the equation of state, consistent with 
all data, is given by 

h(9) = 9- 0.76201(36) 9 s + 8.04(11) x 10" 3 9 5 , (6.3) 

(errors on the last digits in parentheses) that is obtained from p 2 = 2.8667 . This 
expression of h(9) has a zero at 9q = 1.154, which corresponds to the coexistence 
curve. The coefficient of 9 7 in eq. (|67J) is smaller than 10 -3 . Note that for the 
largest value of 9 2 which corresponds to 9 2 , the 9 5 term is still a small correction. 
Finally the corresponding values for the e-expansion are /13 = —0.72, /15 = 0.013. 
These values are reasonably consistent, because a small change in /13 can be 
cancelled to a large extent by a correlated change in p. 



The Widom scaling function /(x), eq. (|3.1|) , can then easily be obtained nu 



merically from h{9) and compared with other determinations. The main dis- 
agreement with other predictions comes from the region x — > 00, i.e. from the 
small magnetization region, where the predictions of the present method should 
be specially reliable. 
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6.3 Amplitude ratios 

Various amplitude ratios can then be derived from h(9) and the values of the 
critical exponents determined earlier from longer series. They involve ratios of 
functions of 9 at 9 = and at 9 = 9q where 9q is the zero of h(9) closest to the 
origin. 

Tables [| and ^| contain a comparison of several amplitude ratios as obtained 
from RG, lattice calculations (High temperature series and Monte-Carlo simu- 
lations) and experiments on binary mixtures, liquid- vapour, uniaxial magnetic 
systems. 



Table | 
Amplitude ratios. 





A+/A- 


c+/c~ 


Rc 


R x 


e — exp., |p7 




0.524 ±0.010 


4.9 




1.67 


e — exp., H 




0.527 ±0.037 


4.73 ±0.16 


0.0569 ± 0.0035 


1.648 ±0.036 


d = 3, g7 




0.541 ±0.014 


4.77 ±0.30 


0.0594 ±0.001 


1.7 


d = 3, || 


1 


0.537 ±0.019 


4.79 ±0.10 


0.0574 ±0.0020 


1.669 ±0.018 


HT series [29 




0.530 ± 0.003 


4.77 ±0.02 


0.0564 ±0.0003 


1.662 ±0.005 


MC [ 31 1 




0.560 ±0.010 


4.75 ±0.03 






bin. mix. 


0.56 ±0.02 


4.3 ±0.3 


0.050 ±0.015 


1.75 ±0.30 


liqu. — vap. 


0.48-0.53 


4.8-5.2 


0.047 ±0.010 


1.69 ±0.14 


magn. syst. 


0.49-0.54 


4.9 ±0.5 







Table | 
Other amplitude ratios. 





Ro 


R 3 


C 4 /C 4 


HT series [5C 

d = 3, [II 
e-expansion, 


), 29 

In 


0.1275 ±0.0003 

0.12584 ±0.00013 
0.127 ±0.002 


6.041 ±0.011 

6.08 ±0.06 
6.07 ±0.19 


-9.1 ±0.2 

-9.1 ±0.6 
-8.6 ± 1.5 



The results from the e-expansion of p7| , |4q| are obtained by direct Pade sum- 
mation of each corresponding series), while the results of are obtained by 
first summing h(9) and then computing ratios, as explained in section 5/2. The 
results from the d = 3 fixed dimension expansion of ref. |I7| refer to direct 
summation up to 0(g 5 ) while the d = 3 values of @ again rely on the method 
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(see also 



explained in ref. fl3[. High temperature results are taken from |p0| , 
[|58| , |59|j ) . Experimental data are extracted from ||60|| . to which we refer for more 
results and references. One notes the overall consistency of the results obtained 
by different methods. 



7 Concluding remarks 



Within the framework of renormalized quantum field theory and renormalization 
group, the presently available series allow, after proper summation, to determine 
precisely critical exponents for the iV-vector model and the complete scaling 
equation of state for 3D Ising-like (N = 1) systems. In the latter example 
additional technical tools, beyond Borel summation methods, are required in 
which the parametric representation plays a central role. From the equation of 
state new estimates of some amplitude ratios have been deduced which seem 
reasonably consistent with all other available data. 

Clearly a similar strategy could be applied to other quantities in a magnetic 
field, in the scaling region. Note also that an extension of the e-expansion of the 
equation of state for N = 1 to order £ 4 or even better e 5 , that does not seem 
an unrealistic goal, would significantly improve the e-expansion estimates and 
would therefore be quite useful. Finally the present approach could be extended 
to systems in the universality class of the {4> 2 )\ field theory for N ^ 1, provided 
expansions of the renormalized effective potential at high enough order can be 
generated. 

More generally it must be emphasized that only a small number of universal 
quantities, as predicted by renormalization group, have yet been calculated. In 
addition to static infinite volume quantities, for which much work remains to 
be done, dynamic properties have not even be touched, finite size calculations 
would be useful for comparison with computer simulations. 

Considering the increase in computer power more perturbative calculations 
will become feasible. What has been demonstrated here is that once the series 
are available, summation methods have been developed which lead to precise 
determinations. 



Acknowledgments. The author thanks R. Guida for careful reading of the 
manuscript. 
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